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Using non-equilibrium Green's function formalism, we investigate nonlinear transport and charg- 
ing effects of gated graphene nanoribbons (GNRs) with even number of zigzag chains. We find a 
negative differential resistance (NDR) over a wide range of gate voltages with on/off ratio ~ 10 6 
for narrow enough ribbons. This NDR originates from the parity selection rule and also prohibition 
of transport between discontinues energy bands. Since the external field is well screened close to 
the contacts, the NDR is robust against the electrostatic potential. However, for voltages higher 
than the NDR threshold, due to charge transfer through the edges of ZGNR, screening is reduced 
such that the external potential can penetrate inside the ribbon giving rise to smaller values of off 
current. Furthermore, on/off ratio of the current depends on the aspect ratio of the length/width 
and also edge impurity. Moreover, on/off ratio displays a power law behavior as a function of ribbon 
length. 
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I. INTRODUCTION 

Graphene is one of the intriguing new materials, which 
has been studied extensively since Novoselov et alP fabri- 
cated it by micromechanical cleavage in 2004. In fact, flat 
structure of graphene makes its fabrication more straight- 
forward than carbon nanotubes. Moreover, dreams of 
carbon nanoelectronic approach to the reality based on 
planar graphene structures. This structure overcomes 
some difficulties of nanoelectronics based on carbon nan- 
otubes, by using lithography, one-dimensional ribbon 
patterns on graphene sheets^. For achieving realistic 
nanoelectronic applications based on graphene nanorib- 
bons (GNR), width of ribbon have to be narrow enough 
that a transport gap is openecP^. Sub-10 nm GNR field- 
effect-transistors with smooth edges have been obtained 
in Ref.[l] and demonstrated to be semiconductors with 
band-gap inversely proportional to the width and on/off 
ratio of current up to 10 6 at room temperature. 

The origin of transport gap which is opened in a gate 
voltage region of suppressed nonlinear conductance is 
still not well understood^!. Based on the tight-binding 
approach, GNRs with armchair shaped edges are ei- 
ther metal or semiconducto^^^. Moreover, in this 
approach, zigzag edge ribbons are metal regardless of 
their widths^. While ab initio calculations^ predict that 
regardless of the shape of the edges, GNRs are semi- 
conductor. Two factors are responsible for transport 
gap: the edge disorder leading to localization^! and the 
confinementF"^. However, in nonlinear regime, trans- 
port gap is also opened by transition selection rules which 
originates from the reflection symmetry^. 

Similar to carbon nanotubes, electronic transition 
through a ZGNRs follows from some selection rules. The 
rotational symmetry of the incoming electron wave func- 
tion with respect to the tube axis is conserved while pass- 
ing through nanotubesS^. Correspondingly, the trans- 
verse reflection symmetry of the incoming and outgo- 
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FIG. 1: Gated Zigzag graphene nanoribbon which is divided 
into three regions: left, right and central region. Dotted rect- 
angular is the unit cell which is used for finding the band 
structure of graphene ribbons. Lower panel shows a field- 
effect transistor structure based on graphene ribbons where 
the gate voltage is applied on the whole system. 



ing wave functions results in the parity conservation in 
ZGNRs with even number of zigzag chain^MJI. As a 
consequence of the even-odd effect, a negative differen- 
tial resistance (NDR) region appears in the I-V charac- 
teristic curve of P-N even ZGNR junctions^. A similar 
NDR behavior has been also reported in P-N nanotubc 
junctional Historically, NDR was first observed in the 
degenerated N-P diode junctions^. Nowadays, NDR has 
been reported in many other molecular deviceiP^Ml 

Motivated by the work done in Ref.pQ, in this pa- 
per, we investigate nonlinear transport through gated 
even ZGNRs by using non-equilibrium Green's function 
(NEGF) approach. The NDR region with on/off ratio 
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of the current up to 10 6 appears in the current-voltage 
characteristic curve. This NDR is induced by transport 
gaps which are opened by two selection rules governing 
electron transition through ZGNRs: (i) the parity conser- 
vation, and ( ii) that allowed transition are between con- 
nected bands^. Based on band structure analyzing, we 
show that transport gap opened by the second selection 
rule is filled for ribbons wider than lOnm. So, sub-lOnm 
ribbons with long enough length provide experimental 
manifestation of the NDR phenomenon in I-V curve of 
GNRs. On the other hand, the gate voltage regulates the 
current flow by shifting the blocked energy regions with 
respect to the Fermi level. Moreover, on/off ratio of the 
current displays a power law behavior as a function of 
ribbon length as M 7 5 . 

Our calculations show that the details of the electro- 
static potential profile along the ribbon can not affect 
the emergence of NDR. The same conclusion has been 
reported by Ref . |20j . but they have not elaborated on 
the physical reason behind this robustness. By follow- 
ing the self-consistent charge and potential profiles at 
different voltages, we demonstrate that at low voltages, 
strong screening of the external potential at contacts re- 
sults in a flat electrostatic potential along the ribbon. 
Subsequently, the e-e interaction at a mean field level, 
does not change the magnitude of I on . However, for volt- 
ages higher than the NDR threshold V on , the transfer of 
charge along the edges, leads to more reduction in I Q f{ 
which improves the switch performance. 

This paper is organized as follows: although the 
formalism has been presented elsewhere^, to be self- 
contained, we briefly introduce the Hamiltonian and 
NEGF formalism in section II. In section III, we dis- 
cuss selection rules governing even ZGNRs. The origin of 
NDR seen in the I-V curve is explained in section IV. We 
demonstrate in section V that the e-e interaction does 
not have a significant effect on the phenomena of NDR 
in the I-V curve. The last section concludes our results. 



II. HAMILTONIAN AND FORMALISM 

Fig. ([I]) shows schematic side view of graphene nanorib- 
bon. In presence of source-drain applied potential, rib- 
bon is divided into three regions; left, right electrodes 
and also central interacting region. Gate voltage is ap- 
plied by means of substrate on the graphene plate. The 
interacting Hamiltonian governing the electron dynamics 
is written in tight-binding approximation. This Hamil- 
tonian is a functional of charge density: 

H{n} = £iOi + [(Xi - x )/L - 0.5]V sd + £\ U l3 5n 3 ])c\c t 

+ T, <lJ> t (4c j +c i cl), 

(1) 

£j shows onsite energy of i th carbon atom and t rep- 
resents the hopping integral between nearest neighbor 
atoms. One n orbital is considered per each site for 
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FIG. 2: transmission (right panel) and band structure of right 
(center panel) and left (left panel) electrodes for ZGNR with 
6 unit cells in length and 4 zigzag chains. Applied bias is 
considered to be a) V s d ~ and b) V s d = 1.0. Gate voltage 
is V g = 0.5V. Here, A and 7 are the energy separation of up- 
per/lower group of bands from the central bands at the Dirac 
point (k x = ±27r/3a) and k x — 0, respectively. Moreover, 8 
is the half- width of the central bands at k x =0. 



graphene system. Without losing any generality, we set 
onsite energies (e^) of all sites equal to zero. All en- 
ergies are in units of tc-C = 2.5eV. Application of a 
gate voltage is achieved by shifting atomic onsite ener- 
gies in all three regions. The applied source-drain po- 
tential, Vgd, and the gate voltage, V g , preserve transverse 
symmetry with respect to the ribbon axis (X direction 
in Fig. (I])). Linear variation of the source-drain voltage 
along the ribbon is the solution of the Laplace equation 
with Dirichletc boundary condition on the contacts. Uij 
is the electrostatic Green's function and Sn,i — rij — n® is 
the change in the self-consistent charge rij from its initial 
equilibrium zero-bias value. This third term is the direct 
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Coulomb interaction created by the bias-induced charges 
at a mean field level which is the solution of Poisson 
equation. The electrostatic Green's function for a distri- 
bution of charges between two parallel conducting planes 
located at x — 0,L which are held at zero potential, 
has the following form: 



U(!>,-&') = 2 J™dkJ {ak) 



s'mh(kz < ) sinh(/c(L — z> )) 
sinh(A;L) 



yJ(x-x') 2 + (y-y') 2 + U H 2 



(2) 



where Uh is the Hubbard parameter whose semi- 
empirical value for carbon^ is about Atc-C- This param- 
eter determines the strength of electron-electron interac- 
tion. This electrostatic Green's function is appropriate 
for the kernel of Ohno-Klopmann modeP^. 

First of all, we must find the self-consistent charge and 
electrostatic potential by using NEGF formalism. The 
retarded and advanced Green's function matrix subjected 
to the central portion of the ribbon is as the follows: 



G r > a (E, n) = [(E ± r/)I - H{n} - 



(3) 



where 77 — > + . "7" is the unit matrix 
retarded self-energies due to scattering by the left/right 
electrodes. These self-energies are independent of the 
charge density. To determine the self-energy, one needs 
to calculate the surface Green's function of semi-infinite 
electrodes g p (E) by using the Lopez-Sancho's method^. 
The escaping rate T of electrons to the electrodes is re- 
lated to the self-energies as T p = i[££ — £*]. Hav- 
ing Green's function, one can find total charge (n = 
' q ) by separate calculations of the equilibrium 



and non-equilibrium charges by using the retarded and 
lesser Green's functions, respectively. It is simply demon- 
strated that in the coherent regime, lesser Green's func- 
tion can be represented by the retarded and advanced 
Green's functions which are determined in Eq.([3]), 
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(4) 

where f p = 1/[1 + cxp( k J^ )\ shows the Fermi function 
of the electrodes, and /iq = = fi^ show electrochemi- 
cal potentials of the left and right electrodes. The initial 
charge n® is calculated by the equilibrium integration in 
zero bias. To obtain the charge, coupled equations of 
[T] and [3| a re self-consistently solved by using Broyden's 
method 29 . The transmission coefficient T(E : V) is de- 
fined in terms of self-consistent Green's functions as: 



T = Tr[GT fl G Q r L ] 



(5) 



The current passing through nanoribbon is calculated 
by the Landauer formula at zero temperature^: 
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(6) 
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FIG. 3: Current- voltage characteristic curve for different gate 
voltages when the ribbon size is as (N, M) = (4, 6) . The 
effect of different parameters such as size effect, electrostatic 
potential and also gate voltage is investigated on I-V curve. 
Hubbard parameter (Uh) is on site Coulomb repulsive. 



III. BAND STRUCTURE AND SELECTION 
RULES 



To understand the transport properties of ZGNRs, 
one needs to study the band structure along with trans- 
mission curves. Band structure of a wide ZGNR pos- 
sesses two well-separated valleys (conically-shaped dis- 
persion curves) at K and K points which can be ef- 
fectively described by 2 + 1 Dirac theory. Edge states 
give rise to fiat bands at charge neutrality point. The 
band structure of ZGNR is calculated in tight-binding 
approximation. The eigenvalues and eigenfunctions of 
the system are extracted by diagonalization of the follow- 
ing HamiltoniarPE H = Hq + To_ + ie lk:rX + T a _ie~ lkocX ] 
where Hq is the tight-binding Hamiltonian for the unit 
cell shown in Fig. ([!]). In the nearest neighbor approxi- 
mation, ro, +1,70,-1 are overlap matrices between a given 
cell ('0') and its right, left neighboring cells, respectively. 
These matrices are 2N x 27V dimensions where N is the 
number of zigzag chains. The eigenvalues E(m,k x ) of 
this Hamiltonian, give the energy spectrum of graphene 
ribbon shown in Fig.([2|. Here, m and k x represent the 
band index and longitudinal wave vector, respectively. 
The eigenfunction can be symmetric (even parity) or 
asymmetric (odd parity) with respect to the mirror sym- 
metry axis^l. As a result of this mirror symmetry, all 
eigenfunctions of each band have either '+' or '-' parity, 
which is indicated in Fig.|2]). If N is even, the lowest 
conduction band has even parity and the highest valance 
band has odd parity. If the bands are labeled by the num- 
ber m from bottom to top as m — 1,2, ..,2 AT, parity of 
m th band is (— l)(" l + 1 ). In ZGNR, we consider longitudi- 
nal source-drain potentials I4d which are invariant under 
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mirror symmetry. Therefore, such V s( j, does not destroy 
the parity of the energy bands, although it changes the 
longitudinal momentum of electron during transport. 

Electronic transition is controlled by two selection 
rules. First, the parity is conserved in tunneling of elec- 
tron through even-ZGNRs. Therefore, at zero source- 
drain voltage, one can expect full transmission which is 
shown in Fig.([2ja). In this case, all bands with the same 
parity are energetically aligned and there is no gap in the 
transmission curve. In Fig.([2ja), the energy of transmis- 
sion curve is shifted by the gate voltage (0.5V). Parity 
of each band is indicated by plus/minus signs. In the 
range of 2A, there is one conducting channel which re- 
sults in unit transmission coefficient. Transmission curve 
a t Vsd = could also be extracted from Fig.Q which 
shows a contour plot of transmission in plane of energy 
and V sd . 

Fig.([2]b) represents band structures of electrodes 
which arc shifted with respect to each other due to the 
source-drain voltage Vsd = 1.0V. The gap in the trans- 
mission curve of Fig.|2]b), AB region, indicates that 
transport between bands of opposite parity is blocked. 

The second selection rule which governs electron trans- 
port through ZGNRs, is that, electron transition is 
allowed between connected bands. Figs.([2ja,b) show 
the band structure classified in three different groups; 
namely, central, upper and lower bands which are indi- 
cated by solid, dashed and dashed-dot-dot lines, respec- 
tively. The common feature of bands in each group is 
that, they are connected at the zone boundary, while 
distinct groups are disconnected. When one considers 
the electron transport, the longitudinal momentum k Xl 
of electrons changes as a result of applied V s( \- The pre- 
cise form of this variation in k x , crucially depends on 
profile of the superimposed longitudinal potential. These 
groups are disconnected from each other from the point 
of longitudinal momentum. Variation of momentum of 
electron k x depends on the shape of superimposed longi- 
tudinal potential. The transport properties for smoothly 
varying V^d, are significantly different from V sc j profiles 
with sharp spatial variations. The electronic transition 
between an eigenstate (mi, k) in the right electrode and 
an eigenstate (ra-i , q) in the left electrode is proportional 
to Fourier transform of longitudinal voltage and structure 
factor^, 

(ip mi (fc) I V sd (x) | V™ 2 (q)) = SV sd (k - q), (7) 

where structure factor of S is equal to [1 + ( — l) p ™i+ p '«2] 
for even ZGNRs and parity of band m is equal to P m — 
(— l)( m + 1 ). Parity selection rule in even ZGNRs origi- 
nates from this structure factor. This parity selection 
rule is mesoscopic analogue of chirality factors governing 
transport of Dirac electrons in planar graphene^S. 

In our calculations, we apply a constant gate voltage 
to the whole system without any longitudinal variation. 
As a result, the gate voltage does not change momen- 
tum of electron. However, linear variation of the applied 
source-drain bias (with the slope V s d/L) changes the elec- 
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FIG. 4: Contour plot of transmission with respect to the en- 
ergy and V E d for the system size (N,M) = (4,6) and the applied 
gate voltage V g = 0.5V . Darked Oblique lines shows the cur- 
rent integration window. Points marked by A, B, C, D corre- 
spond to the horizontal lines with the same name in Fig. {[2J . 
Lines of a, /3 are the trace of points with similar names in 
Fig.|3|. Fermi energy was fixed at zero Ef = 0. 

tron momentum. So, smooth variation of the potential 
in longer ribbons results in a small momentum variation 
of electron. Consequently, transition of electron between 
disconnected bands is forbidden when the length of rib- 
bon is so large that one can assume V s d(k — q) — > 6(k — q). 
Therefore, a smooth potential in the longitudinal direc- 
tion can just scatter the electron among the class of states 
belonging to the same group of connected energy bands. 

Now let us focus on the two transport gap regions: AB 
and CD in Fig.([2jb). The AB gap is a consequence of par- 
ity selection rules, while the CD gap is due to blockage of 
transition between disconnected groups. As can be seen 
in Fig.([2jb), the AB gap is proportional to the source- 
drain voltage, V s d- Moreover, this gap is independent of 
the ribbon width. Of course, in wide ribbons, upper and 
lower band groups approach to the central group, espe- 
cially at point k x — 0, where 7 in Fig.([2ja) tends to zero 
as log-normal. 

When the ribbon width is increased, the separation 7 
between the upper/lower and central groups of bands, is 
reduced, which tends to loosen the second selection rule 
based on band groups; hence filling in the gaps. However 
when we increase the ribbon length, our classification 
of bands into connected groups is recovered. Therefore 
the AB gap is essentially governed by the aspect ratio of 
ZGNR. 

The CD gap is equal to A — 5 + V s d, where A and 5 are 
the energy separation of upper/lower group from the cen- 
tral bands at the Dirac point, and the half width of the 
central bands at k x = 0, respectively. The dependence 
of A on width N is: A cx (2.13 ± o.02)iV"( - 864±0 003 ), 
while 5 has a Log-Normal behavior which asymptotically 
approaches to the constant value of 0.9738 ± 0.0002 as 
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FIG. 5: Contour plot of the current in terms of V s d and V g 
for a zigzag graphene nanoribbon with (N,M) = (4,6). Vertical 
dotted line corresponds to V g = 0.5V. a and /3 are those 
points which was shown in Figs( 3|4 1, 



FIG. 6: a) On-off ratio of the current increases as a power 
law with the ribbon length for TV = 4 and V g = 0.51/. b) 
On-off ratio of the current decreases with the ribbon width 
(circle and triangular points) while it is disappeared for edge 
impurity (diamond points) stronger than 0.3tc-c (inset fig- 



TV goes up to 10. The conducting region BC in Fig.([2j 
can exist only, when 5 — A < V s d < 5. From the de- 
pendence of A and 5 on N, the CD gap exists if N is 
less than 30. Hence, NDR is estimated to be observable 
for ribbon width < 7nm. The lowest achieved ribbon 
width is sub-lOnm-wide (~ 2 ± 0.5nm]^ with the length 
~ 236nm. Such long ribbons with small width provide 
fascinating experimental manifestation of the selection 
rules in transport properties. 



IV. NEGATIVE DIFFERENTIAL RESISTANCE 

Fig.Q shows current- voltage characteristic curve of a 
ZGNR with 4 zigzag chains and 6 unit cells in length. In 
the case of zero gate voltage, flow of current is blocked 
due to the parity selection rule, while at a given V s d, 
gate bias turns the current on. After a range of V^d in 
which the current remains unchanged, current begins to 
reduce with increasing V s d- NDR threshold voltage V on 
decreases with gate voltage for V g < 0.6V. Dependence 
of NDR threshold voltage on the gate voltage can also be 
seen in Fig. ^ . This NDR also symmetrically appears in 
the negative polarity of Vsd ■ The NDR threshold voltage 
and ion remain unchanged in the presence of the electron- 
electron interaction (with a given Hubbard term U = 
Atc-c)- However, reduction of the current in off state, 
Tog, is intensified when one takes electrostatic potential 
into account. 

To understand the origin of NDR, it is helpful to look 
at the 3D contour-plot of transmission in plane of energy 
and V s d which is presented in Fig.Q. Blocked energy 
intervals AB and CD, which are indicated in Fig.Q, cor- 
respond to those intervals shown in Fig. For voltages 



lower than the vertical line a (V < V a ), transmission is 
a nonzero constant for the whole region of the conduc- 
tion window represented in Eq.Q. As a result, current 
increases proportionally to V s d- In the voltage interval 
[Va, Vp], the blocked region AB, originating from the par- 
ity selection rule, contributes to the current integration 
window of Eq.([6|. However, nonzero range of transmis- 
sion remains unchanged along with V s d resulting in the 
fixed current in the voltage range [V a , Vp\. So, current re- 
mains unchanged in this range. For voltages V > Vp, the 
CD gap contributes in the current integration window, 
and consequently the NDR phenomenon emerges. 

Regarding the importance of gate voltage in the cur- 
rent flow, let us investigate the effect of the gate volt- 
age on / — V s d curve by contour plotting of the current 
with respect to V s< i and V g in Fig.Q. For gate volt- 
ages \V g \ < 0.1V, shift of transmission is not remarkable 
enough to contribute to conducting channels in the cur- 
rent integration. So, current is blocked by the parity 
selection rule. In the range 0.1V < \V g \ < 0.6V, contri- 
bution of conducting region BC in transport is accom- 
panied with the blockage arising from AB and CD gaps 
in voltages V > V on . As a consequence, current reduces 
after a threshold voltage. In this range, on/off ratio of 
the current increases and V on reduces with increasing the 
gate voltage. 

As can be seen in I-V curves of Fig. ([3]), off current 
reduces for longer ribbons which enables us to achieve 
high performance switches by increasing aspect ratio of 
the ribbon length/width. The reason is connected to 
smooth variation of the applied potential along the rib- 
bon such that during transport, electrons are scattered 
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among those states which belong to continuous bands. 
As a consequence, blockage originating from electronic 
transition between disconnected bands is intensified by 
increasing the ribbon length. In fact, when the length of 
ribbon increases, transmission in the AB and CD gaps 
decreases exponentially. 

Since off current is induced by contribution of the gaps 
in the current integration, I s efficiently decreases with 
increasing the ribbon length (M). Fig.(|6]a) shows that 
lon/IoS displays a power law behavior as a function of 
the ribbon length for large M: I on /I s oc M v where 
r] = 7.5061 ± 0.03505. As an example, for M = 50, 
on/off ratio goes up to 10 6 which suggests experimental 
fabrication of high performance switches based on the 
GNR nanoelectronics. 

Experimentally, it was observed in Ref. [3] that the 
room-temperature on/off ratio induced by the gate volt- 
age increases exponentially as the GNR width decreases. 
They observed that I on /I s is equal to I, 5, 100 and > 10 s 
for W = 50nm, 20nm, lOnm and sub-lOnm, respectively. 
Similarly, as shown in Fig.(6]b), on/off ratio calculated 
for the set up considered in this paper, also decreases with 
the ribbon width, while reduction of on/off ratio can be 
compensated by considering longer ribbons. However, 
NDR phenomenon is disappeared for the ribbons wider 
than 7nm. 

In ab initio calculations 6 , by using hydrogen- 
termination of zigzag edges, mirror-symmetry of ZGNRs 
and consequently parity conservation could be retained. 
Correspondingly, by several repetition of the heat treat- 
ments and hydrogenation, it is also possible to create 
well-ordered H-terminated edges in experiment^. How- 
ever, the edge states with energies about —0.1 to 0.2 
eV have been experimentally observed^ that emerge at 
hydrogen-terminated zigzag edges. To simulate the edge 
states and the effect of symmetry breaking on NDR phe- 
nomenon, it is assumed to dope one of the ZGNR edges 
by slight impurity. Edge impurity is considered to ap- 
ply as a change in the on-site energy of the edge atoms 
(e Q ) with respect to on-site energy of the other atoms. 
In case of edge disorder, e a plays the role of the averaged 
on-site energy of the edge atoms. Inset figure indicated 
in Fig.([6jb) shows that on/off switching reduces with the 
edge impurity strength, however, NDR still emerges for 
e a < 0.3t C -c- 



ELECTROSTATIC POTENTIAL AND 
CHARGING EFFECT ON NDR 



Emerging phenomenon of negative differential resis- 
tance in I-V curve is not destroyed by the e-e interaction 
and is independent of the details of electrostatic poten- 
tial profile. However, interaction reduces off-current as 
shown in the I-V curves of Fig. To substantiate the 
above claim, comparison of transmission curves in the 
presence and absence of the e-e interaction is useful. It is 
apparent from Fig.([7ja) that for voltages less than V^n, 
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FIG. 7: Transmission curve of zigzag graphene nanoribbon 
with 4 zigzag chains in width for voltages a) 0.75 volt b) 1.25 
volt. Interacting and non-interacting curves are compared 
with each other. 



transmission in conducting channels is robust against the 
e-e interaction while transmission increases in the gaps 
with respect to the non-interacting case. But this en- 
hancement is slight enough which can not affect the emer- 
gence of NDR. However, for voltages V > V on , interaction 
lowers transmission coefficient in the conducting channels 
(such as BC region) in which higher subbands participate 
in transport. Such behavior is corroborated in Fig.([7|b) 
which indicates transmission at that voltage correspond- 
ing to the off-current, V ff . Reduction in the transmission 
coefficient of the conducting channels results in further 
reduction of the off-current. To explain the reason for 
such phenomenon, it is necessary to study the poten- 
tial and charge profiles. Electrostatic potential averaged 
on each unit cell is represented in Fig.([8ja) in terms of 
the ribbon length. For voltages less than V on , potential 
sharply drops only at the contact regions which connects 
the system to electrodes. In such a case, external poten- 
tial is strongly screened by redistributed electrons and, 
electrostatic potential of the central atoms remains close 
to zero. Screening is performed by discharging of elec- 
trons from the area connected to the source and their ac- 
cumulation around the drain electrode. These facts are 
obvious from transferred charge and electrostatic poten- 
tial profiles represented in Figj9j Since U(n — n ) deter- 
mines electrostatic behavior of the potential, discharging 
of electrons weakens the external potential penetrated 
from the source electrode. Moreover, charge accumula- 
tion around the drain electrode prevents potential drop 
in the central part of the system. So in the case of strong 
screening, potential drops only at the contacts. However, 
when the applied bias goes beyond 1 volt, screening is be- 
ing weakened and external potential can penetrate inside 
the central region. The reason as to why screening is 



7 



(a) 



(b) 



«J 0.25- 



O 
0. 



p 



£ -0.25 » 



O 
O 



LU 



edge charge 
middle charge 
total transferred charge 




Unit Cell Position( A) 



0.5 1 1.5 2 
Source-Drain Voltage(V sd ) 



1 



0.14 
0.12 
0.11 
0.10 
0.08 
0.07 
0.06 
0.04 
0.03 
0.01 
0.00 
-0.01 
-0.03 
-0.04 
-0.05 



0.63 
0.54 
0.45 
0.36 
0.27 
0.18 
0.09 
0.00 
-0.09 
-0.18 
-0.27 
-0.36 
-0.45 
-0.54 
-0.63 



Transferred Charge Profile for 1 .25 Volt 



• • • • • 

Potential Profile for 1 .25 Volt 



FIG. 8: a) Edge, middle and total transferred charge in terms 
of source-drain applied bias, b) Electrostatic potential per 
unit cell in terms of unit cell position for source-drain voltages 
0.5, 0.75, 1.25 volts. The gate voltage is for all curves equal 
to V a = 0.5 V. 



FIG. 9: Transferred charge and Electrostatic potential profiles 
for V a d = 1.25 volt and V g = 0.5 volt in the weak screening 
case. Due to charge transfer through the ribbon edges, screen- 
ing is so weak that the external potential can penetrate inside 
the central portion. 



weak, can be sought in charge distribution. Figure (|8jb) 
illustrates that for voltages less than 0.5 volt, in and out 
flow of charge are balanced with each other such that the 
total transferred charge remains close to zero. However, 
around the voltage V on and voltages above, the charge 
is mainly transferred from the edges of the ribbon, so 
that the source electrode does not inject further charge 
to middle atoms of the ribbon. As a consequence, by 
increasing the applied bias and so gradient of the poten- 
tial along the central region, charge depletion is mainly 
enhanced in the middle bar area of ZGNR. On the other 
hand, since the only way for transporting electrons is the 
edge atoms, significant accumulation of charge appears 
along the two edge lines of ZGNR. 

In summary, at voltages less than V on , electrostatic po- 
tential is only dropped at the contacts and therefore mo- 
mentum of electrons is only varying in the area where the 
potential drops, while longitudinal momentum of elec- 
tron remains unchanged across the central portion. In 
other words, potential steeply drops in the low-area dis- 
trict around the contacts which results in violation of 
the blockage rule which governs on transition between 
disconnected energy bands. Subsequently, transmission 
coefficient slightly increases in the blocked energy ranges. 
In other words, in this case, an increase in gradient of the 
potential facilitates electronic transport in the blocked 
energies. Note that interaction preserves transverse sym- 
metry, so the parity selection rule still governs electronic 
transport. Therefore, the AB gap induced by the parity 
conservation still survives for voltages larger than V^n. 
For voltages V > V on , electrostatic potential gradually 
penetrates into the whole system so that the potential of 
the central region is not flat. In addition, because the 



edge transport of electrons dominates, the transverse po- 
tential is deeper in the middle of ZGNR than its edges. 
Therefore, the band structure of the interacting central 
region differs from the band structure of electrodes. As a 
consequence, for voltages V > V on , transmission of con- 
ducting channels and also off current reduces. 



VI. CONCLUSION 

In this paper, using non-equilibrium Green's function 
approach, we investigated nonlinear transport and charg- 
ing effects of graphene nanoribbons with even zigzag 
chains. Current flow is controlled by the gate voltage 
applied on the whole of sub-lOnm ribbons. In this range 
of widths, two selection rules govern on the electronic 
transition are: (i) the parity conservation, and (ii) al- 
lowed transition between connected bands. As a result, 
a negative differential resistance (NDR) in TV curve is 
appeared in the presence of the gate voltage. Further- 
more, on/off ratio of the current increases with the rib- 
bon length as a power law behavior up to 10 6 , while rib- 
bon width and edge impurity reduce on/off ratio. Emer- 
gence of the NDR phenomenon is not sensitive to details 
of the electrostatic potential profile. Because of strong 
screening in low voltages, the major potential drop takes 
place at the contacts. However, in voltages larger than 
the NDR threshold, due to charge transfer through the 
ribbon edges, screening is so weak that the external po- 
tential can penetrate inside the central portion. As a 
consequence, off current reduces in comparison to non- 
interacting ribbons. 



VII. ACKNOWLEDGES and Dr. S. A. Jafari for a critical proof-reading of the 

manuscript and for his useful comments and suggestions. 

We wish to acknowledge Prof. K. Esfarjani for 
collaboration^ in the development of NEGF formalism 



* Electronic address: cheraghchi@dubs.ac.ir 

1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

2 Z. Liu, K. Suenaga, P. J. F. Harris, and S. Iijima, Phys. 
Rev. Lett. 102, 015501 (2009), B. ozyilmaz, P. Jarillo- 
Herrero, D. Efetov, D. A. Abanin, L. S. Levitov, and P. 
Kim, Phys. Rev. Lett. 99, 166804 (2007). 

3 M. Y. Han, B. ozyilmaz, Y. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007); F. Molitor, A. Jacobsen, C. 
Stampfer, J. Gttinger, T. Ihn, and K. Ensslin, Phys. Rev. 
B. 79, 075426 (2009). 

4 X. Li, X. Wang, L. Zhang, S. Lee, H. Dai, Science, 319, 
1229 (2008); X. Wang, Y. Ouyang, X. Li, H. Wang, J. Guo, 
H. Dai, Phys. Rev. Lett. 100, 206803 (2008). 

5 F. Molitor, A. Jacobsen, C. Stampfer, J. Gttinger, T. Ihn, 
and K. Ensslin, Phys. Rev. B. 79, 075426 (2009). 

6 Y. W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 
97, 216803 (2006). 

7 F. Sols, F. Guinea, and A. H. Castro Neto, Phys. Rev. 
Lett. 99, 166803 (2007). 

8 K. Nakada, M. Fujita, G. Dresselhaus, and M. S. Dressel- 
haus, Phys. Rev. B. 54, 17954 (1996). 

9 L. Brey, H. Fertig, Phys. Rev. B. 73, 195408 (2006); Phys. 
Rev. B. 73, 235411 (2006). 

10 H. Zheng, Z. F. Wang, T. Luo, Q. W. Shi, J. Chen, Phys. 
Rev. B. 75, 165414 (2007). 

11 L. Malysheva and A. Onipko, Phys. Rev. Lett. 100, 186806 
(2008). 

12 E. R. Mucciolo, A. H. Castro Neto, and C. H. Lewenkopf, 
Phys. Rev. B. 79, 075407 (2009). 

13 Z. Li, H. Qian, J. Wu, B-L. Gu, and W. Duan, Phys. Rev. 
Lett. 100, 206802 (2008). 

14 A. A. Farajian, K. Esfarjani and Y. Kawazoe, Phys. Rev. 
Lett. 82, 5084 (1999). 

15 A. Cresti, G. Grosso, G. P. Parravicini, Phys. Rev. B. 78, 
115433 (2008); Phys. Rev. B. 77, 233402 (2008). 



16 A. R. Akhmerov, J. H. Bardarson, A. Rycerz, and C. W. 
J. Beenakker, Phys. Rev. B. 77, 205416 (2008). 

17 Z. Li, H. Qian, J. Wu, B. L. Gu, W. Duan, Phys. Rev. 
Lett. 100, 206802 (2008). 

18 J. Nakabayashi, D. Yamamoto, and S. Kurihara, Phys. 
Rev. Lett. 102, 066803 (2009). 

19 K. Wakabayashi and T. Aoki, Int. J. Mod. Phys. B. 16,4897 
(2002). 

20 Z. F. Wang, Qunxiang Li, Q. W. Shi, Xiaoping Wang, 
Jinlong Yang, J. G. Hou, Jie Chen, Appl. Phys. Lett. 92, 
133114(2008). 

21 L. Esaki, Phys. Rev. 109, 603 (1958) 

22 D. Dragoman, M. Dragoman, Appl. Phys. Lett. 90, 143111 

(2007) . 

23 H. Cheraghchi, K. Esfarjani, Phys. Rev. B. 78, 085123 

(2008) , see also its references. 

24 S. E. Ulloa, G. Kirczenow, Phys. Rev. B. 35, 795 (1987). 

25 J. D. Jackson, Classical Electrodynamics, 2 nd Ed , pl32. 

26 K. Esfarjani and Y. Kawazoe, J. Phys.: Cond. Matt. 10, 
8257 (1998). 

27 K. Ohno, Theor. Chim. Acta 2, 219 (1964); G. Klopman, 
J. Am. Chem. Soc. 86, 4550 (1964). 

28 M. B. Nardelli, Phys. Rev. B. 60, 7828 (1999). 

29 C. G. Broyden, Math. Comput. 19, 577 (1965), also See, 
e.g., K. Ohno, K. Esfarjani, and Y. Kawazoe, Computa- 
tional Materials Science from Ab Initio to Monte Carlo 
Methods (Springer, Berlin, 1999). 

30 G. C. Liang, et al, Phys. Rev. B. 69, 115302, sec IIA 
(2004). 

31 M. Ezawa, Phys. Rev. B. 73, 045432 (2006). 

32 D. Rainis, F. Taddei, F. Dolcini, M. Polini, R. Fazio, Phys. 
Rev. B. 79, 115131 (2009). 

33 Y. kobayashi, K.-I, Fukui, T. Enoki, and K. Kusakabe, 
Phys. Rev. B. 73, 125415 (2006); Y. Niimi, T. Matsui, H. 
Kambara, K. Tagami, and H. Fukuyama, Phys. Rev. B. 
73, 085421 (2006). 



